Renormalization group for evolving networks 
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We propose a renormalization group treatment of stochastically growing networks. As an exam- 
ple, we study percolation on growing scale-free networks in the framework of a real-space renormal- 
ization group approach. As a result, we find that the critical behavior of percolation on the growing 
networks differs from that in uncorrelated nets. 
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Evolving networks with a complex distribution of con- 
nections attract much interest from a wide circle of re- 
searchers [1-5]. The physicists' contribution to the rapid 
progress in this field is based on understanding that net- 
works are objects of classical statistical mechanics and 
can be effectively studied by using standard approaches 
of statistical physics. 

In this communication, we demonstrate how a real- 
space renormalization group approach (or positional 
renormalization group), which is traditional in critical 
phenomena theory (see, e.g., Refs. [6-11] and references 
therein), can be applied to stochastically growing net- 
works. For the demonstration, we consider a bond per- 
colation problem, although other models of cooperative 
behavior on evolving networks can be investigated in a 
similar way. 

Percolation on uncorrelated equilibrium networks is 
well studied [12-15]. In fact, many other cooperative 
models (with discrete symmetry of the order parameter 
and absence of frustrations) on equilibrium networks (the 
spread of diseases [16,17], the Ising and the Potts models 
[18], etc.) show a behavior similar to that for the per- 
colation. In this respect, percolation problems are very 
representative. 

Percolation on a growing network is a more complex 
problem than for equilibrium ones. The reason for this 
complexity is a wide spectrum of correlations, which are 
inevitably present in growing nets. Note that the corre- 
lations between the degrees of the nearest-neighbor ver- 
tices [19-21], which were thoroughly studied in equilib- 
rium networks [22-25], is only a particular kind of the 
correlations. 

We study the following problem. First we grow an infi- 
nite size network and then consider a classical bond per- 
colation problem on it. That is, randomly chosen edges of 
the infinite network are simultaneously removed. A frac- 
tion p of edges is retained. We use the model of a stochas- 
tically growing, undirected, highly-clustered, scale-free 
network [5] (see Fig. 1), which is ideally suited for a 
real-space renormalization group procedure. In principle, 
the percolation problem for this network can be exactly 
solved, at least, at some particular case. We, however, 
use this model for the first demonstration of the renor- 
malization group method for growing random networks. 



The second moment of the degree distribution P(k) of 
this growing network diverges (degree is the total number 
of connections of a vertex, sometimes it is called "con- 
nectivity"). We show that the "percolation threshold" 
is zero, p c — 0, that is the percolating cluster (the gi- 
ant connected component) exists at each p ^ (super- 
stability against random damage). The 7 exponent of 
the degree distribution P(k) oc k 1 of the net is below 3. 
We find that all the derivatives of the relative size of the 
giant component M(p) over p diverge at the percolation 
threshold. Moreover, M(p < 1) ~ e - const / p . T h is differs 
from the corresponding critical singularity for percola- 
tion on uncorrelated scale-free networks with exponent 
7 < 3. 

The network is constructed in the following way [5]. 
The growth starts from a single edge connecting two ver- 
tices (t = 0). At each time step, each edge of the network 
transforms as shown in Fig. 1: 

(a) with probability q, an edge "creates" a vertex which 
is attached to both the end vertices of its "mother" 
edge, or 

(b) with the complementary probability 1 — q, this edge 
creates a bare vertex. 





FIG. 1. Edge transformation which generates the stochas- 
tically growing scale-free network. At each time step, each 
edge of the network transforms into one of the two shown 
configurations with the complementary probabilities q and 
1-q. 
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Of course, we can equally leave an edge unchanged in 
item (b), and, at first sight, the creation of bare vertices 
seems superfluous. However, this inessential feature of 
the transformation will be convenient for us. 

In the particular case q = 1, this network turns out 
to be the simple deterministic graph that was introduced 
in Ref. [4] and studied in detail in Ref. [26]. This graph 
(q = 1) belongs to a wide class of deterministic grow- 
ing nets which are currently being intensively studied 
[27-29]. If q < 1, the growth is stochastic. As q — > 
(but q 0) the network coincides with a network grow- 
ing by attaching a new vertex to the ends of a randomly 
chosen edge at each time step [30] . The latter has practi- 
cally the same degree distribution as the Barabasi- Albert 
model (7 = 3) but is strongly clustered. Indeed, the net- 
work that is generated by the transformation shown in 
Fig. 1 has numerous loops of length three. So it has a 
large clustering coefficient at any value of q. 

One can easily modify the transformation in such a way 
that the resulting network will turn out to be exactly the 
Barabasi- Albert model in the limiting case q — > (q ^ 0) . 
Let at each time step, each edge creates a new vertex 
which, with probability q, is attached to one end vertex 
of the edge and, with the same probability, is attached 
to the other end vertex. At q = 1, we again obtain the 
deterministic graph [4,26], and as q — > 0, the network ap- 
proaches the Barabasi- Albert model with zero clustering. 

We, however, apply the transformation that provides 
the strongly clustered networks at each q. Note that we 
use transformations generating networks with the small- 
world effect [1], that is the average shortest-path length 
of the network grows logarithmically with the network 
size. 

Simple calculations show that the degree distribution 
of the network is scale-free with exponent 



7 = 1- 



ln(l + 2g) 
ln(l + q) 



(1) 



The spectrum of degrees is continuous at q < 1. As q de- 
creases from 1 to 0, 7 grows from 1 + In 3/ In 2 = 2.585 . . . 
to 3. The average number of edges in the network grows 
as (1 + 2qY, where t is the number of a time step. 

For demonstration, we use the simplest, rather naive 
version of the real-space renormalization group approach 
[6-9,11,32]. Let us outline the procedure in application 
to our network. We stop the grows at a time step T — > 00 
and spoil the network: if an edge is present in the undam- 
aged net, then this edge is retained in the damaged net 
with the probability p. Then we invert the transforma- 
tion in Fig. 1 and define n = T — t for the inverted 
transformation, which is actually a decimation proce- 
dure. Further, we introduce the probability p n that if 
an edge connects a pair of vertices of the undamaged net 
at t = T — n, then at the n-th step of the decimation for 
the damaged net, there is connection between these ver- 
tices. Here, pa = p. One can easily derive the following 
recursion relation for p n : 



Pn+l = q[Vn + 3^(1 - p n ) + Pn(l ~ Pn) 2 } 

+ (1 - q)Pn ■ (2) 

Its structure is evident from Fig. 1. 

Let us find the dependence of the relative size M of 
the percolating cluster on p. Here, for convenience, we 
define M as the fraction of edges that belong to the gi- 
ant connected component of the network (the percolating 
cluster). So, M(p = 1) = 1 in the undamaged network. 
As is usual, for calculating M, we find the average num- 
ber m n of edges in the "percolation" configurations of the 
renormalization group transformation at an n-th step: 

P n m n = 

q[3 ■ pl_ x + 2 • 3p£_!(l - Pn-l) + 1 • Pn-i(l - Pn-i) 2 } 

+ (1 - q)pn-l (3) 

[compare Eqs. (2) and (3); notice also the multiple p n 
on the left-hand side of Eq. (3)]. Then, M(p) can be 
obtained from the following relation 



M 



m. 



=1 1 + 2 <1 ' 



(4) 



where we have taken into account that the number of 
edges in the undamaged network increases by \ + 2q times 
at each time step of the evolution. 

The final forms of the recursion relation for p n and the 
expression for the relative size of the percolating cluster 
are 



Pn+l = Pn[qPn{t ~ Pn) + 1] , 

1 + 2qp n (2-Pn) 

n=l 



M(p) =pH 



l + 2q 



(5) 



where po = p. One can see that the recursion relation for 
p n has only two fixed points, and 1. At any p = p ^ 0, 
p n approaches 1 as n — > oo, which indicates the presence 
of the percolating cluster. So, the percolating cluster of 
this network cannot be eliminated by the random removal 
of edges at any value of the parameter q. 

By using the relations (5), we numerically obtain the 
dependence M(p) (see Fig. 2). From Fig. 2(a), one can 
see that the curves M(p) weakly depend on the param- 
eter q if p is large enough. The essential difference is 
visible in the range of small p, see Fig. 2(b). The analyt- 
ical analysis of the relations (5) immediately yields the 
asymptotic behavior of M{p) at small p: 



M(p «1)~ exp 



ln(l + 2g) ! 



(6) 



We do not write down the preexponential factor, since 
our simplified, demonstrative renormalization group pro- 
cedure certainly cannot give its proper form. In fact, at 
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each step of the procedure, the complex spectrum of the 
probabilities of the realization of various edge configura- 
tions is renormalizcd to (one can also say, is substituted 
by) a single delta-function distribution. If a cooperative 
model under consideration does not contain frustrations, 
then usually, this substitution is not dangerous (we dis- 
cuss only qualitative results). This is the case for per- 
colation. Only for q -C 1, this approximation can cause 
serious problems [e.g., very inaccurate values of the re- 
sulting factor in the exponential of Eq. (6)]. 
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FIG. 2. (a) Relative size of the percolating cluster, M, as 
the function of the fraction of retained edges, p. The curves 
correspond to the values of the parameter q: 1.0, 0.8, 0.6, 0.4, 
0.2, 0.1, 0.05, from top to bottom, (b) The same as (a) but 
in the region of small p. 



There are two possible reasons for this difference. The 
first is the wide spectrum of correlations induced by the 
growth. The second is the high clustering and numerous 
loops in our network. The loops, in principle, may lead 
to fluctuation effects. (These effects are absent in coop- 
erative models on uncorrelated networks, which have a 
tree-like local structure.) 

One has to admit that a real-space renormalization 
group approach has a reputation of an uncontrolled ap- 
proximation. Moreover, in our demonstration, we have 
used the simplest version of the approach, which has al- 
lowed us to obtain analytical results. In principle, one 
can use more refined versions of the real-space renormal- 
ization group procedure (e.g., see Ref. [10]), where, how- 
ever, numerical calculations are necessary [33]. 

We considered one of traditional cooperative models 
on undirected growing networks. However, the renor- 
malization group approach can also be used as a tool 
for studying structural properties (e.g., a distribution of 
loops [34], etc.) of evolving networks. Directed networks 
can be considered in a similar way as for randomly di- 
rected percolation on a lattice [11]. Furthermore, the 
transformation that we used for generating the growing 
random network (see Fig. 1) is only a simple example 
and can be easily generalized [5] . 

In this communication, the renormalization group ap- 
proach was applied to evolving networks. One should 
mention the applications of a real-space renormalization 
group technique to static networks: to small-world net- 
works in Ref. [35] and, recently, to a random network of 
masses connected by springs in Ref. [36]. 

In summary, we have applied the real-space renor- 
malization group procedure to stochastically growing 
networks. In the framework of this approach, we have 
considered the percolation problem for scale-free grow- 
ing networks with exponent 7 in the interesting region 
7 < 3. The percolation threshold is zero, p c — 0, but the 
observed critical behavior essentially differs from that 
for uncorrelated networks. Our results demonstrate new 
possibilities of a renormalization group method. 
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Expression (6) shows that all the derivatives of M(p) 
diverge at the "percolation threshold" p c = 0. The type 
of this divergence differs from that for uncorrelated scale- 
free networks with 2 < 7 < 3, where M(p <C 1) ~ 
pV(3-7) anc l only if 7 = 3, the relative size of the perco- 
lating cluster behaves as M(p <1)~ e -const/ P ^ gee r^ e f 
[15]). 
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